(submitted to Phys. Rev. Lett.) 

First-principles statistical mechanics study of the stability of a sub-nanometer thin 
surface oxide in reactive environments: CO oxidation at Pd(100) 
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We employ a multiscale modeling approach to study the surface structure and composition of 
a Pd(100) model catalyst in reactive environments. Under gas phase conditions representative of 
technological CO oxidation (~ 1 atm, 300-600 K) we find the system on the verge of either stabilizing 
sub-nanometer thin oxide structures or CO adlayers at the surface. Under steady-state operation 
this suggests the presence or continuous formation and reduction of oxidic patches at the surface, 
which could be key to understand the observable catalytic function. 

PACS numbers: 82.65.+r, 68.43.Bc, 68.43.De, 68.35.Md 



Under realistic operating conditions heterogeneous cat- 
alysts are exposed to reactant pressures of the order of at- 
mospheres. In recent years, there has been an increasing 
awareness that for oxidation catalysis, i.e. under atmo- 
spheric O2 conditions, the surface of the employed tran- 
sition metal (TM) catalysts might get oxidized. Instead 
of the pristine TM surface often studied in atomic detail 
under ultra- high vacuum (UHV) conditions, it could then 
be the surface oxide that only forms under reactive con- 
ditions which actuates the catalytic activity under atmo- 
spheric steady-state operation. Such situation induced 
by the reactive environment has recently been revealed 
for CO oxidation at so-called Ru catalysts, with two note- 
worthy conclusions 0,0,0: (i) There is a material's gap, 
i.e., in the initial period when the Ru metal is exposed to 
the reactive environment an oxide film is formed at the 
surface; and (ii) there is a pressure gap, i.e., on the once 
formed Ru02 surface CO and O are only adsorbed in the 
proper ratio (necessary for efficient catalysis), if pressure 
and temperature of the reactants are sufficiently high. 

As an element from the middle of the TM series Ru 
forms rather strong chemical bonds and Ru02 is in fact 
the stable bulk phase under the reactive conditions 
This is different for the late transition and noble met- 
als, which are equally or even more frequently used in 
oxidation catalysis. In this study we therefore set out 
to further elucidate the situation by using CO oxida- 
tion over the Pd(100) model catalyst as an example. 
For this system, in-situ reactor scanning tunneling mi- 
croscopy (STM) experiments had interpreted an observed 
roughening of the surface as oxide formation, accompa- 
nied with a notable increase in the catalytic activity p. 
The first-principles statistical mechanics calculations de- 
scribed in this Letter fully confirm this interpretation, in 
that indeed oxide formation in the reactive environment 
plays a key role also at the more noble Pd model cata- 
lyst. The striking gist is that this oxide is not a bulk-like 
film that once formed is stable and actuates the cataly- 
sis. Instead, our investigation points at the relevance of 
a sub-nanometer thin surface oxide film, which is likely 



continuously formed and reacted away in the on-going 
reaction. 

Aiming at quantitative atomic-scale insight we employ 
a multiscale modeling approach to analyze the surface 
structure and composition of Pd(100), when exposed to a 
reactive gas phase consisting of O2 and CO. The atomic- 
level energetics of the system is obtained from density- 
functional theory (DFT) calculations, performed within 
the highly accurate full-potential augmented plane wave 
+ local orbit als (LAPW/APW+lo) scheme 6]. The ex- 
tended surface is described within a slab geometry and 
we carefully checked that the numerical uncertainty in- 
troduced by the finite basis set, as well as by the su- 
percell parameters does not affect the reported results 
[3 . The only notable approximation with respect to the 
energetics underlying our approach is thus the approxi- 
mate treatment of the exchange-correlation interaction, 
for which we use the generalized gradient approximation 
(GGA-PBE) [3, and we will critically discuss this point 
below. 

In a first step, we account for the effect of a finite 
gas phase environment within a constrained atomistic 
thermodynamics approach 2]. Here, the DFT energet- 
ics enters into the calculation of the Gibbs free energy of 
adsorption of possible states of the system. Neglecting 
kinetic effects of the on-going catalytic reactions on the 
surface structure and stoichiometry (at this stage), the 
surface is considered to be in equilibrium with i sepa- 
rate reservoirs representing the i gas phase species, each 
characterized by partial pressures pi and temperature T 
(or more conveniently by a chemical potential A/i ? ;(T,p^) 
summarizing this two-dimensional dependence 2]). We 
verified that for the focus of this work, bulk and sur- 
face vibrational and configurational contributions to the 
Gibbs free energy of adsorption are not significant Q and 
correspondingly neglect them here. 

Within this approach we are then in a position to com- 
pare the stability of a wide range of possible surface struc- 
tural models under (T,po 2 ,Pco)-conditions ranging from 
UHV to atmospheric pressures and up to elevated tern- 
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FIG. 1: (Color online) Surface phase diagram for the Pd(100) surface in "constrained thermodynamic equilibrium" with an 
environment consisting of O2 and CO. The atomic structures underlying the various stable (co-) adsorption phases on Pd(100) 
and the surface oxide, as well as a thick bulk- like oxide film (indicated by the bulk unit-cell), are also shown (Pd = large spheres, 
O = small spheres, C = white spheres). Phases involving the bulk oxide are cross-hatched, phases involving the surface oxide 
are hatched. The dependence on the chemical potentials of O2 and CO in the gas phase is translated into pressure scales at 
300 K and 600 K. The thick black line marks gas phase conditions representative of technological CO oxidation catalysis, i.e. 
partial pressures of 1 atm and temperatures between 300-600 K. The three vertical grey (green) lines correspond to the gas 
phase conditions employed in the kinetic Monte Carlo simulations shown in Fig. El (see text). 



peratures. Aiming to assess a possible oxide formation 
at the surface, we consider bulk PdO to represent thick 
bulk-like oxide films |^ . In addition we exploit the recent 
characterization of a (y/E x y/E)R27° (henceforth called 
y/E) surface oxide structure, which corresponds to a sub- 
nanometer thin film of PdO(lOl) on the surface In- 
cluding a large pool of combinatorially representative or- 
dered (co-) adsorption phases of O and CO on this y/E 
surface oxide and on the pristine Pd(100) surface, we ar- 
rive at a set of 191 possible surface configurations. Out 
of these, we find only 11 to be most stable in a certain 
range of O2 and CO chemical potentials, or equivalently 
(T,po 2 5.Pco)-gas phase conditions. Figure [I] summarizes 
these stability regions in (A/iQ, A/ico)-space, where we 
additionally include the corresponding pressure scales for 
two representative temperatures. 

A natural starting point to discuss the overall topology 
of this surface phase diagram is the bottom left part of 
Fig. This corresponds to vanishing concentrations of 
gas phase species and correspondingly, the clean Pd(100) 
surface results as the most stable system state. Moving 
from this situation to the right in the surface phase di- 
agram, e.g. by increasing the O2 pressure while keeping 
the CO pressure low, we find surface structures contain- 
ing an increasing amount of oxygen. These are a p(2 x 2) 
adsorption phase with O atoms adsorbed at on-surface 
hollow sites at a coverage © = 0.25 monolayer (ML), the 



y/E surface oxide (O = 0.8 ML) and bulk-like thick oxide 
films; a sequence that was recently confirmed in in-situ 
surface X-ray diffraction measurements [llj . In the com- 
plementary case of low oxygen pressure and increasing 
CO pressure (moving from the bottom left in Fig. Q 
upwards) a series of ordered CO adsorption phases on 
Pd(100) are stabilized, most of which have already been 
characterized in UHV experiments For simultane- 

ously increased O2 and CO gas concentrations, i.e. mov- 
ing now up right along the diagonal in Fig. d there are 
three phases with O or CO adsorbed at the y/E surface 
oxide structure, but interestingly no co-adsorption struc- 
tures on Pd(100). Although our set of surface configu- 
rations compared in the thermodynamic approach com- 
prises also a large number of such structures, we find 
them largely disfavored by overall repulsive lateral inter- 
actions among the adsorbates, which is consistent with 
the experimentally observed tendency of adsorbed O and 
CO to form separate domains on the Pd(100) surface [l^ . 

The overall topology of the constrained surface phase 
diagram can thus be divided into three different parts. 
First, a region of gas phase conditions, where bulk-like 
thick oxide films are stable (cross-hatched area). Sec- 
ond, a region comprising various adsorption phases on 
the y/E surface oxide (hatched area), and finally, a region 
with different adsorption phases on Pd(100). Particu- 
larly intriguing is that gas phase conditions representa- 
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tive of technological CO oxidation catalysis (pi ~ 1 atm, 
T ~ 300 — 600 K) fall almost exactly on the boundary 
between the stability region of adsorption phases on the 
y/E surface oxide and the stability region of CO-covered 
Pd(100) phases, cf. Fig. E] On the basis of these results 
we would correspondingly be able to draw two important 
conclusions with respect to the relevance of oxide forma- 
tion in the reactive environment: (i) The formation of 
a thick bulk-like oxide film at the surface as in the case 
of the Ru catalyst can be ruled out for Pd(100) under 
any gas phase conditions representative of technological 
CO oxidation catalysis, (ii) The stability region of the re- 
cently characterized y/E surface oxide, on the other hand, 
does extend to such conditions. 

Before further discussing the implications of these find- 
ings, it is appropriate to scrutinize the major approxima- 
tion of our approach, namely the xc functional employed 
in the DFT calculations. To this end, we recomputed Fig. 
[fusing the local-density (LDA and another gradi- 
ent corrected (GGA-rPBE [0|) xc functional, which both 
yield notably different O2 and CO gas phase binding en- 
ergies compared to the GGA-PBE 15]. In both cases we 
obtain the same topology of the surface phase diagram, 
i.e. we find the same stable phases in the same regions 
as in Fig. Q However, in general there are notable shifts 
in the location of the boundaries between the various 
phases in (A/io, A/ico)-space. Strikingly, this does not 
apply to the boundary between the CO-covered Pd(100) 
phases and the adsorption phases on the y/E surface ox- 
ide. For all three xc functionals, the catalytically rele- 
vant gas phase conditions lie very close to this boundary, 
but still within the stability region of the surface oxide. 
This situation is neither changed, when approximately 
accounting for the bulk and surface vibrational contribu- 
tions to the Gibbs free energy of adsorption 0, Q . Con- 
sidering the configurational entropy contributions would 
smear out the infinitely sharp boundary drawn in Fig. 
H over a small, but finite range of chemical potentials. 
While the width of this phase coexistence range can be 
estimated Q , we defer its discussion to the kinetic simu- 
lations below, where such configurational effects are fully 
included. 

The constrained thermodynamics employed up to here 
is understood as a starting point, helping us to iden- 
tify the regions in (T,pi) -space that deserve closer at- 
tention by a higher-level approach considering chemi- 
cal reactions and kinetics. In a second step, we there- 
fore explicitly account for such effects by performing 
first-principles kinetic Monte Carlo (kMC) simulations 
to specifically investigate the stability range of the y/E 
surface oxide. In these simulations we follow the time 
evolution of the system as described by a master equa- 
tion and use rates for the individual elementary processes 
that are obtained from the DFT GGA-PBE energetics to- 
gether with transition-state theory |3| . We focus here on 
the lattice structure of the y/E surface oxide, in particular 
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FIG. 2: Kinetic Monte Carlo (kMC) results for the stabil- 
ity of the y/E surface oxide. Shown is the average occupation 
O hoi of the hollow sites with oxygen in the y/E lattice struc- 
ture. The chosen gas phase conditions (T = 300, 400, 600 K, 
po 2 = 1 atm, variable CO pressure) correspond to the three 
vertical lines drawn in Fig. [I] In the intact y/E surface oxide 
structure G Q hoi = 100 %, and the reduction obtained in the 
kMC simulations occurs at CO partial pressures quite close 
to those corresponding to the stability boundary as obtained 
in the constrained thermodynamic approach (transition from 
the hatched to the plain areas in Fig. 0, marked here by the 
corresponding vertical lines. 



on the surface kinetics involving the hollow and bridge 
sites. We consider all non-concerted adsorption, desorp- 
tion, diffusion and Langmuir-Hinshelwood reaction pro- 
cesses involving these sites. Dissociative O2 adsorption 
requires two neighboring sites, and is only hindered by 
a computed dissociation barrier of 1.9 eV in the case 
of adsorption into two bridge sites. CO adsorption is 
unimolecular and not hindered by adsorption barriers. 
Nearest-neighbor lateral interactions are taken into ac- 
count in the elementary process rates and we verified 
that the kMC set-up reproduces the part of the thermo- 
dynamic phase diagram involving all y/E phases (hatched 
part of Fig. HJ, when the reaction events are not allowed 
to occur p. 

The kMC simulations were performed for fixed 
(T,po 2 ?Pco)-conditions on a lattice comprising (50 x 50) 
surface unit-cells, evaluating the average occupation at 
the various surface sites after steady-state was reached. 
We carried out simulations in the temperature range 300- 
600 K and for fixed po 2 = 1 atm. With an initially set 
low CO partial pressure, we are then in the middle of 
the thermodynamic stability region of the y/E surface ox- 
ide shown in Fig. [l] And indeed, we exactly obtained 
average occupations pertinent to the geometric structure 
of the corresponding phases in Fig. namely two on- 
surface oxygen atoms sitting in their preferred hollow 
sites 10], while the bridge sites are either empty or filled 
with one or two CO molecules or one O atom per sur- 
face unit-cell. Under such gas phase conditions, kinetic 
effects due to the on-going reactions are thus negligible 
and we simply reproduce the results obtained within the 
atomistic thermodynamic approach. 

In subsequent simulations we then gradually increased 
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the CO partial pressure, which brings us closer and closer 
to the boundary of the stability region of the y/E surface 
oxide as indicated by the vertical lines in Fig. ^ The 
most important surface reactions that can occur between 
the adsorbed species under these gas phase conditions are 
O hoiiow + C0 bridge _^ C q 2 with a calculated reaction bar- 
rier of A^ e h a c 1+cobr = 0.9 eV and O brid § e + CO hollow 
CO2 with A£^Qb r c +cohol = 0.5 eV. Kinetic effects when 
these reactions consume surface species faster than they 
can be replenished from the gas phase may lead to the 
destabilization of the y/E surface oxide at already lower 
pco than predicted by the stability boundary in the ther- 
modynamic phase diagram. To determine the onset of 
this structural decomposition in the kMC simulations we 
focus on the average occupation O hoi of O species at 
the hollow sites. As just discussed, these sites are always 
fully occupied in the intact surface oxide structure and 
we interpret a Oohoi reduced to 80% as a sign for the on- 
set of the structural decomposition. The corresponding 
data as a function of the CO partial pressure for various 
temperatures is shown in Fig. 

The depopulation of the hollow sites occurs at all tem- 
peratures over a certain range of CO partial pressures. 
The deduced critical pressures for the onset of oxide 
decomposition depend nevertheless little on the present 
choice of reduction to 80 % occupation, and we also ver- 
ified they would be little affected by slight 0.1 eV uncer- 
tainties in the employed DFT GGA-PBE surface reaction 
barriers. As apparent from Fig. [2 the stability range of 
the y/E surface oxide defined by these critical pressures is 
in all cases quite similar to the one obtained within the 
constrained atomistic thermodynamic approach. Still, 
the small variation of the stability region induced by the 
reaction kinetics has a notable effect. Whereas in the 
constrained equilibrium approach the catalytically rele- 
vant environments (pi = 1 atm, 300-600 K, black line in 
Fig. ^) were always fully within the stability region of 
the y/E surface oxide, this is no longer the case when the 
reaction kinetics are taken into account. At the lower 
temperatures considered, slightly oxygen-rich conditions 
are now required to stabilize the surface oxide film. At 
400 K, this is a partial pressure ratio of po 2 /pco ~ 10 : 1, 
which is in good agreement with the recent in-situ reac- 
tor STM experiments p performed at such gas phase 
conditions. At the higher temperatures, the stability of 
the surface oxide is enhanced, and our simulation results 
predict that at 600 K, the surface oxide is still stable even 
up to a pressure ratio of 1:1. This suggests that at least 
at these higher temperatures the catalytically most rele- 
vant gas phase conditions correspond to a transient sit- 
uation where either sub-nanometer thin oxide structures 
or CO adlayers are just stabilized at the surface, such 
that small changes (or even fluctuations) in the partial 
pressures may have profound consequences on the global 
or local surface composition and structure. 

Employing DFT together with constrained thermody- 



namics as well as statistical mechanics we therefore pro- 
vide evidence for the relevance of oxide formation at the 
Pd(100) model catalyst under gas phase conditions rep- 
resentative of technological CO oxidation. Depending on 
temperature, the system is on the verge of stabilizing 
sub-nanometer thin oxide structures or CO adlayers at 
the surface. Similar to the case of CO oxidation at Ru 
catalysts, it is therefore again a fluctuation-rich region 
close to boundaries in the computed surface phase dia- 
gram that is of particular interest for the catalytic func- 
tion. While in Ru the coexistence is between different 
adsorption phases on an otherwise stable oxide film, at 
the more noble Pd(100) it involves the presence of CO 
on pristine Pd and CO on surface oxide patches. Under 
steady-state operation this suggests possible oscillations 
between these two states, which could then also drive a 
spatio-temporal pattern formation. 
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